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Liverpool, L69 3BX, United Kingdom 

We review the application of the critical point large Nf self-consistency method to 
QCD. In particular we derive the 0(1/Nf) d-dimensional critical exponents whose 
e-expansion determines the perturbative coefficients in MS of the field dimensions, 
/3-function and various twist-2 operators which occur in the operator product ex- 
pansion of deep inelastic scattering. 



1 Introduction 

The renormalization group equation, (RGE), plays an important role in com- 
paring predictions made in a quantum field theory with observations of na- 
ture. The fundamental ingredients in the RGE are the renormalization group 
functions. Since these are rarely known exactly even for the simplest of field 
theories one has to be content with approximate perturbative solutions; the 
accuracy being dependent upon how many orders in the perturbative coupling 
constant one can compute the RGE functions. This is a highly technical and 
tedious exercise partly because the number of Feynman diagrams at even one 
loop can sometimes be excessive. Also the results depend on how one removes 
the ultra-violet infinities. For theories which particle physicists are interested 
in such as quantum chromodynamics, (QCD), which is the gauge theory de- 
scribing the strong interactions, roost high order calculations of these functions 
are performed in the MS schemeErta For instance, the /3-function of QCD has 
been deduced at third order in this scheme. Recently information on various 
scattering, amplitudes has been produced at the same level in an impressive set 
of papersu Due to the complexity of such calculations, having independent and 
alternative methods to check the high order structure of the RGE functions is 
important. 

One such method has been made available through the properties of the 
RGE in the neighbourhood of a fixed point which is defined to be a non- 
trivial zero of the /3-function. There it is known that the critical exponents 
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which characterize the phase transition correspond to the functions of the RGE 
evaluated at the critical coupling. So if one can compute exponents directly 
then information on the RGE functions is obtainable^ This has been achieved 
in impressive articles by Vasil'ev et al for the 0(N) a modelQ There critical 
exponents were determined in arbitrary dimensions order by order in powers of 
1 /N when N is large. Those results are in total agreement with the e-expansion 
at the fixed point of the same exponents deduced explicitly at 5-loops in MS. 
Garnered by that success it is therefore a worthwhile exercise to develop the 
1/Nf method for QCD, where Nf is the number of quark flavours, in relation 
to the present state of the art calculations. 



2 Basic ideas 



We recall the basic ideas for deducing arbitrary dimensional critical exponents 
in the 1/A^jexpansion. First from the two loop /3-function of QCD in d- 
dimensions[±l there is a fixed point at 



3e 



T(R)N f AT 2 (R)N 



f 
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where d = 4 — 2e. If, for example, a general RGE function takes the form 

7(ff) = c l9 + (c 2 N f + d^g 2 + (c 3 N 2 + d 2 N f + ei )g 3 + 0(g 4 ) (2) 

where the coefficients {q, cij, &j . . .} are independent of Nf, then the associated 
exponent at leading order in 1/Nf is 

-. 00 

^9c) = w 5>[3e/7W + 0(l/N 2 ) (3) 

f r=l 

So provided "f(g c ) can be computed directly in the large Nf limit its e-expansion 
gives the leading order sequence of coefficients {ci} of 7(3). 

The exponents are defined with reference to the action of the theory one 
is interested in. For QCD this takes the form 

L = ~ (4) 

where is the quark field, is the gluon field, = + T a A a ^ G^ v = 
d^Al - d v A a + f abc A b Al, Tfj is the generator of the colour group whose 
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structure constants are f abc , I < i < Nf, 1 < I < N c and 1 < a < (N 2 - 1). 
The canonical dimensions of the fields of Eq. 4 at g c are defined by demanding 
that the action is dimensionless. The anomalous dimensions are defined to be 
the extra portion of the full dimension of the field or operator and essentially 
are a measure of the effect of radiative corrections. For instance, in the pealing 
region where the propagators of Eq. 4 behave in the limit k 2 — > oo, asp 



rinv - (1 - b) 



(5) 



where A and B are momentum independent amplitudes and b is the covariant 
gauge parameter, we define 

a = A f - 1 + \n > P = 1 - V ~ X (6) 

with d = 2/i. Here \ IS the dimension of the quark gluon vertex operator 
and rj is the quark anomalous dimension. Expressions for these anomalous 
dimensions are deduced from studying the scaling dimensions of the next tg 
leading order corrections to the 2 and 3 point Green's functioruusing Eq. 50 
For an arbitrary gauge parameter, the leading order results arep 

= C 2 (fl)[(2 M -l)( M -2)+//%? 

V (2 M -l)( M -2)T(ii)iV / 

, = C 2 (C?)[(2 M -l) + fe( / /-l)]7 ? ° 

V X 2(2/x- l)(n-2)T(R)N f [ ' 

where 77° = - (2/i - 1)(2 - M )r(2 M )/[4r 2 (/i)r(/i + l)r(2 - /*)]. 

In computing these results, which agree with 3-loop perturbative calcu- 
lations in the Landau gaugeja we made use of another well known feature of 
critical point theory. Ordinarily more than one model can be used to de- 
duce exponents at a fixed point and such models are said to be in the same 
universality class. A well known example is the equivalence of the O(N) a 
model and-O(N) 4> 4 theory in three dimensions. For the present case QCD is 
equivalents at leading order in 1/Nf to a non-abelian version of the Thirring 
model, (NATM), which is renormalizable in strictly two dimensions. Its la- 
grangian is 

(A a ) 2 

L = - (9) 

where A is the coupling constant which is dimensionless in 2 dimensions. Elim- 
inating the auxiliary spin-1 field A° yields a 4-fermi term. The benefit of using 
this model, Eq. 9, is that it has a simpler structure to Eq. 4 as the 3 and 4 
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point gluon self interactions are absent. So one need only consider diagrams 
built with the quark gluon interaction. It was shown, thoughp that in the 1/Nf 
limit the 4-fermi model correctly reproduced the 3 and 4 point gluon Feynman 
rules in the approach to four dimensions. In other words with Eq. 9 the effect 
of the 3-point gluon interaction is contained in the graphs with a quark loop. 
This feature occurs implicitly in the calculations we report on later. Further 
in using a covariant gauge, ghost fields have to be included in each lagrangian 
but they give no contribution at leading order. 

3 /^-function 

With this basic formalism the 0(1/ Nf) correction to the QCD /3-function can 
be computedlilJ Ordinarily this is the first step in determining 0(1/ Nj) infor- 
mation as it will encode the next order correction to g c to all orders in e. To 
determine this we compute the related exponent ui = — /3'(g c )/2. It is deduced 
from the last term of Eq. 4 which gives the scaling law 

w = rj + x + XG (10) 

where \G is the critical dimension of the composite operator G — (G^„) 2 when 
computed as an insertion in a Green's function in the non-abelian Thirring 
model. For QED u was originally deduced in 1/Nf by explicithz_performing the 
MS renormalization with an infinite chain of electron bubblesLil The extension 
to the non-abelian case is simpler in the critical approach. Three 2-loop and 
one 3-loop graphs need to be evaluated which are illustrated in Fig. 1. The 




Figure 1: Graphs for 0(X/Nt) contribution to w. 

first two graphs correspond to the QED sector, whilst the remaining two would 
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be absent by Furry 's theorem in QED as their colour group factor is C2(G). 
Consequently, using the critical propagators we find 



{ft -2) - [(2 / x-3)(/i-3)(7 2 (i?) 

uv7„r<m r,9 

(ii) 



(4/x 4 - 18^ 3 + 44^ 2 - 45^+ 14) G 2 (G)] r# 



4(2^-1)^-1) 



The e-expansion of Eq. ll_correctly reproduces the 0(1/Nf) coefficients of 
the 3-loop MS /MunctionEHj With this agreement we can deduce several new 
higher order coefficients. Using the notation 



/2 11 \ 

(3(g) = (d-4)g + (~T(R)N f - -C 2 (G) g 2 + ^ a r [T(R)Nf] 

ra (12) 

for the large Nf leading order part of the /3-function, thcnliJ 
a 4 = - [154G 2 (i?) + 53G 2 (G)]/3888 

a 5 = [(288C(3) + 2U)C 2 (R) + (480C(3) - 229)G 2 (G)]/31104 
a 6 = [(864C(4) - 1056C(3) + 502)G 2 (i?) 

+ (1440C(4) - 1264C(3) - 453)G 2 (G)]/233280 
a 7 = [(3456C(5) - 3168C(4) - 2464C(3) + 1206)G 2 (i?) 

+ (5760C(5) - 3792C(4) - 848C(3) - 885)G 2 (G)]/1679616 (13) 



4 Twist-2 operators 

With the impressive progress that has been made at 3-loops in MS in the 
renormalization of the twist-2 operators of the operator product expansion 
used to understand processes in deep inelastic scatteringcl it is important to 
have some large Nf results available for comparison. Similar to the /3-function 
calculation the critical exponents corresponding to the anomalous dimensions 
of such operators are deduced by inserting the operator into a Green's function 
in the NATM and determining the scaling behaviour of the integrals. The 
operators which we consider are, 

_ i«-is^yi£)/M .. . D^Tfjip J -trace terms 
C^i-Mn _ -n-i^/yn^ ... I}^ 7 - trace terms 

•n-2 

= -—StxG a ^ v D^ . . . D' in - 1 G a J?" - trace terms 
where S denotes symmetrization on the Lorentz indices. 
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For the fermionic twist-2 flavour nonsinglet and singlet operators, ns and 
O s , we deduce at leading order in 1/Nf respectively!! 2 ] 



( „) _ 2C 2 (i?)( M -l) 2 7 ? » 



{2vL-l)(n-2)T(R)N f 



( „) _ (ti-l)C 2 (R) V ? 

Vs 



(2vL-l)(n-2)T(R)Nf 



(n-l)(2/i + n-2) 2 M ' 

(M + n-l)(/x + n-2) + W 
2( M -l)(n-l)(2 M + n-2) f 



(/i + n- l)(// + n-2) 
- ^r(n - l)[(n 2 +n+ 2/i - 2) 2 + 2(/x - 2)(n(n - l)(2/z - 3 + 2n) 
+ 2(/x - 1 + n))]r(2/z)/[(M + n - l)(/i + n - 2)r(2/i - 1 + n)] (14) 

where n is the operator moment, ^>(n) — ip(p — 1 + n) — an d tp( x ) is 

the logarithmic derivative of the T-function. One feature of the singlet sector 
is that the operators do not mix since the gluonic and fermionic operators 
have different canonical dimensions at g c . By contrast in the perturbative 
calculation there is mixing and one has to compute a matrix of anomalous 
dimensions. To compare the e expansion of Eq. 14 with perturbative results 
one realises that in the large Nf calculation the result contained in Eq. 14 is 
in fact the anomalous dimension of the predominantly fermionic eigenoperator 
of the perturbative mixing matrix. Therefore if one computes the eigenvalues 
of the mixing matrixaEj and evaluates them at g c the coefficients of both e 
expansions ought to be in agreement. We record this occurs exactly at the 
3-loop level at leading order in 

More explicitly we present the n-dependence of the coefficient C3, in the 
notation of Eq. 2, of both the nonsinglet and singlet leading order large Nf 
part of the anomalous dimensions. Having the explicit dependence is important 
since the inverse Mellin transform of the anomalous dimensions with respect to 
n determine the Altarelli Parisi splitting functions. These are a function of the 
conjugate variable, x, which is the momentum fraction carried by the partons 
contained in the nucleons, and are in effect a measure of the probability tha)_a 
parton fragments into other partons. First, we have for the nonsinglet casepJ 

of = ls 3 (n) - i p 2{ n) - l Sl (n) + g 



[12n 4 + 2n 3 - 12n 2 - 2n + 3] 
27n 3 (n + l) 3 



(15) 



where Si(n) = 1/r' . To compare with- the results of the explicit 3- 

loop MS calculation for the first few momentspl we have evaluated Eq. 15 for 
various n and presented them in Table 1. For the singlet sector we can deduce 
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n 


„ns 
c 3 


2 


28 
243 


4 


384277 
1944000 


6 


80347571 


333396000 


8 


38920977797 

1 Z ± / ±VZ (XJ { zUUU 


10 


27995901056887 
95850016416000 


12 


65155853387858071 


210582486065952000 


14 


68167166257767019 


210582486065952000 


16 


5559466349834573157251 
16553468064672354816000 


18 


19664013779117250232266617 
56770118727793840841472000 


20 


6730392290450520870012467 
18923372909264613613824000 


22 


16759806821032136669044226177 


46048135637404510767879321600 



Table 1: Values of c ns for various n. 



the n-dependence of the 3-loop coefficient ofj-the anomalous dimension of the 
predominantly fermionic eigenoperator. It isp3 

4 = §*.(») - f 7 S 2{ n) - + g 

2{n 2 +n + 2f[S 2 (n) + Sl(n)} 
3n 2 (n + 2)(n+l) 2 (n- 1) 

- 2S l {n)[l6n 7 + 74n 6 + 181n 5 + 266n 4 + 269n 3 + 230n 2 

+ 44n - 24]/[9(n + 2) 2 (n + l) 3 (n - l)n 3 ] 

- [100n 10 + 682n 9 + 2079ti 8 + 3377n 7 + 3389n 6 

+ 3545n 5 + 3130n 4 + 118n 3 - 940n 2 - 72n 

+ 144]/[27(n + 2) 3 (n + l) 4 n 4 (n - 1)] (16) 
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Similar to 03 s we have evaluated Eq. 16 for low moments and presented the 
results in Table 2. These are in exact agreement with the first four moments of 
the explicit three loop MS results after diagonalizing the mixing matrix Q and 
extracting the leading order large Nf piece corresponding to the dimension of 
the predominantly fcrmionic cigenoperator. 



n 


„s 
c 3 


2 





4 


121259 
720000 


6 
8 


3166907 
13891500 

1328467729 
5038848000 


10 


304337312935261 
1054350180576000 


12 


842357166098254633 


2737572318857376000 


14 


42512567719680559 
131614053791220000 


16 


755896148277147625515451 
2251271656795440254976000 


18 


1121815282809553973842772849 
3235896767484248927963904000 


20 


78640886458671664340562623 
220772683941420492161280000 


22 


4248342909129791924572989157741 
11650178316263341224273468364800 



Table 2: Values of C3 for various n. 

Aside from agreeing with explicit perturbative results up to three loops, 
there are several other checks on the exponents arising from general principles. 
First, as the operators are physical their anomalous dimensions are gauge inde- 
pendent. We have therefore computed Eq. 14 with a non-zero covariant gauge 
parameter b and observed its cancellation in assemblying the contributions 
from the relevant Feynman diagrams in each exponent. Second, for certain 
values of n the corresponding operators reduce to conserved physical currents. 
Provided the conservation of these currents is not spoiled by an anomaly then 
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their anomalous dimensions must be zero to all orders in perturbation theory. 
For the nonsinglet sector the n = 1 case relates to charge conservation, whilst 
the singlet operator with n = 2 corresponds to the energy momentum tensor. 
Therefore for both these respective values the critical exponents of Eq. 14 
must vanish. It is an easy exercise to verify this. Indeed the zero entry for n 
= 2 in Table 2 is a reflection of this general result in the three loop case. 



5 Conclusions 

The critical renormalization group ideasB have proved useful in giving some 
insight into the structure of the MS perturbative coefficients at higher orders 
in QCD. Although we have concentrated on the four dimensional theory the 
results have been expressed as functions of d. Therefore we can also obtain 
information on the three dimensional model. For example, from Eq. 11 

1 10C 2 (G) l\ 

Higher order 1/Nf calculations are possible too. For instance, in the abelian 
sector the dimension of the mass operator, iptfi, is available in <i-dimensions_ 
So when d — 3 the gauge independent electron mass anomalous dimension ispJ 



32 64[3tt 2 - 28] 



~fm(9c) = - — - -^KZI^ + O I TJW 1 (18) 




Such results will be useful for comparing with numerical results for the same 
quantity computed by other methods. Indeed exponents which are known to 
similar orders in other models like the O(N) 4-fermi model and evaluated for 
low N have been in good agreement with lattice resultsLLa 
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